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xo > moreover, the relative and absolute errors converge to as x — > 
oo. These conclusions are based on corresponding results for the digamma 
function i/) := F'/r. Relations with other relevant results are discussed, as 
well as the corresponding computational aspects. This work was motivated 
by studies of exact bounds involving the Student probability distribution. 
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1. Summary and discussion 

My interest to the Wallis ratio 

r(a;+l) 



W{. 



r(x + i) 



was stimulated by recent work [12, 11] on exact bounds involving Student's 
probability distribution, with the density function fp defined by the formula 

rCP±l'] / T,2\-(p+i)/2 
:= 1 + ^ 

for all real x; here p is a positive parameter referred to as the number of degrees 
of freedom. Let us extend this definition by continuity to p = od, so that foo is 
the probability density function of the standard normal distribution: fooix) = 
g-x /2 Jqj. g^jj j.g^j Thus, 

m . Mw(l^) ,L1, 



/oo(0) V P 

for all p € (0, cxj); equivalently. 



Wix) ^r{2x + l)^x+^ (1.2) 

for all X £ (—5, 00). Therefore, bounding the Wallis ratio W{x) is equivalent to 
bounding the ratio r{p), which is what we shall do, as the bounds on r{p) are 
slightly simpler to express and easier to operate with. The ratio r(jj) may also 
be slightly more natural, since r(p) — >■ 1 as p — > 00. 

Theorem 1.1. Take any p e (0,oo). Then one has the identities 

r[p) = UUp)=LUp), (1-3) 



whe 



C/fc(p) :=exp^ 2-i-"Vp,™, (1.4) 



m— 1 



^(-l)^ (^. ) ln(p + j), (1.5) 



and k G {0, 1, . • . , 00}; the series J2m=i ™crp.,„ converges uniformly in p G 
\po,oo), for any po G (0,oo). (^As usual, the sum of an empty family is defined 
to be 0; so, Uo{p) = 1-) 
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Moreover, for all k G {0, 1, . . . } the bracketing inequalities 
Lk{p) < Lk+iip) < r{p) < Uk+i{p) < Ukip) 
hold. For the relative errors 

Uk{p)-r(jp) Uk{p)-r{p) 



(1.7) 



Sup,k{p) ■■= 
SioMp) ■= 



r{p) 



r{p) 



and 



r{p)-Lk{p) Lk{p)~r{p) 



r{p) 



r(j)) 



of the upper and lower approximations one has 

< Sio.kip) < S^pAp) < cxp(p*^fc) - 1, 

where 



1 



Pp,k ■ — 



2k+i 



(1.8) 



(1.9) 



P^ ' ' P 

which converges to as k ^ oo (no slower than geometrically) uniformly in 
p S [po, oo), for any pq G (0, oo), and also converges, for each k G {0, 1, . . . }, to 
as p ^ oo. 

The necessary proofs are deferred to Section 2. 

Note that the definition (1.6) of the lower bound Lk{p) originates in the 
identity 

r{p)r{p + l) = ,/-^ (1.10) 



(1.11) 



\l p+l 

for all p > 0, which easily follows from r(a; + 1) = xT{x). Equivalently, 



W{x)W{x- 



for all X' > — i. Thus, given an upper bound on W{x) (or on r(p)) , one automat- 
ically has a corresponding lower bound — and vice versa; this is a well-known 
trick. 

An important feature of the lower and upper bounds Lk{p) and Uk{p) is that 
they are all algebraic. Therefore, as is done in [11], one can use the well-known 
result by Tarski [17, 4] to solve, in a completely algorithmic manner, systems of 
equalities/inequalities involving such bounds. Let us list here the initial 5 upper 
bounds on r(p); 



Usip) = 



(l-h2a)^ 



\1/16 /(i 



(l + a)ii(l + 3a) 



(l-t-2a)i6(l+4a)\i/32 



-ha)26(l-h3a)e 



where a stands for I /p. Recall that the corresponding lower bounds can be 
obtained by (1.6). 
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Recalling also (1.2) and assuming the correspondence 

(-i,oo) 9 2^ =x < — >p = 2a; + l G (0,cx3), (1.12) 

introduce 

Uk{x)^Uk{p)^Wi and Cu{x) = Lk{p)^Wi. (1.13) 

Then one immediately has the following corollary of Theorem 1.1: 

Corollary 1.2. Relations (1.3) and (1.7) can he rewritten as 

Ckix) < Ck+iix) < C^{x) =. Wix) =U^ix) <Uk+iix) <Uk{x) (1.14) 

for all X E (— i,oo) and k £ {0,1,...}. Moreover, the relative errors of the 
bounds lAk (x) and Ck (x) converge geometrically and uniformly on any interval 
of the form [a;o,oo) for any xq > — ^; also, the relative errors converge to as 
x — > oo. Furthermore, it follows from (1.8)-(1.9) and the inequality W{x) < 

Uq{x) ~ \J X + \ that, for each fc G {0, 1, . . . }, the absolute errors Uk{x) — W{x) 
and W{x) — Ck{x) converge to 0, too, as x ^ oo. 

Theorem 1.1 is based on the following proposition, which provides geomet- 
rically convergent sequences of rational upper and lower bounds on 'ip{x + 1) — 
'(/'(a; + ^) and thus may be of independent interest. Here, as usual, tjj denotes 
the digamma function, that is, the logarithmic derivative of F: 



Proposition 1.3. Take any p E (0, oo) and k E {0, 1, . . . }. Then 

lnr\ 
dp 



d In rin) 

lk{p) < ik+lip) < iooip) = = Uoo{p) < Uk+lip) < Uk{p), (1.16) 



where, for k € {0, 1, . . . , oo}, 



r{y + m) 

J) = 



and 

m— 1 

y'-'--U(y^^^ r(,) 

is the Pochhammer symbol, with y^^^ = 1, corresponding to the general conven- 
tion that the product of an empty family is defined as 1 . Equivalently, 

ik{x) < ik+lix) < ioo{x) = -ipix + 1) - 1p{x + |) = Uooix) < Uk+l{x) < Uk{x) 

(1.18) 

for all X E (— |, oo) and k G {0, 1, . . . }, where, for all k £ {0, 1, . . . , oo}, 
ik{x) := i + 24(p), and Ukix) := i + 2uk{p), 
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again assuming the correspondence (1.12). Moreover, one has the following 
hounds on the errors: 

, , dlnr(p) dlnr(p) „ , , 2-'=-i(A; + 1)! 2-^-^ , 

which converges to as k oo (no slower than geometrically) uniformly in 
p € [poi oo), for any po S (0, oo), and also converges, for each k S {0, 1, . . . }, to 
as p — > oo . 

Let us now go back to the bounds Uk{p) and Li;{p), defined in (1.4) and (1.6). 
It is not hard to see that they follow some simple recursive relations: 

Proposition 1.4. Take any p > and k G {1,2,...}. Then 



"M=^U„^^i,)L„-^iA (1.20) 



where 

Vk{p) -.^ ^Uk{p) and Mk{p):=^Lk{p). (1.22) 

The "initial" conditions are Uq{p) ~ 1, Lq{p) = \J ^o{p) = v^' '^'^'^ 
M,{p) = 

One possible way to apply Proposition 1.4 is to use the first recursive relation 
in (1.21) to compute Vk{p), and then get Uk{p) and Lk{p) by (1.22) and (1.6). 
In view of (1.4) and (1.5), one can also rewrite Uk{p) as 

Uk{p) = {YliP + jY-'^'"^^")'^' ^ , (1.23) 

where Hj^k ■= J2m=ivj 2'^^"'(™) can be found recursively: 



Hj^j = I{j ^ 1} and H^^k+i = ^H^^^ 



k+l 
3 



for nonnegative integers j and k such that j ^ k; here, as usual, I{^} denotes 
the indicator of an assertion A, so that 1{A} equals 1 if ^ is true and equals 
otherwise. Recalling also the known recursive relation (2.1) for the binomial 
coefficients, one sees that the calculation of the bound Uk{p) takes 0{k'^) oper- 
ations such as addition and multiplication, followed by raising the result to the 
power 1/2^^+^; here we used the fact that is takes 0(log7V) multiplications to 
raise a given real number to the power of a natural number N . Since the relative 
errors 5\o,k{p) and 5up,k{p) are 0(1/2*^), it takes 0(log^ operations to make 
the relative errors less than eg (0, 1). 
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Moreover, observe that the integers Hj^k do not depend on p. Therefore, 
in rather typical situations (such as plotting or numerically integrating) when 
one would need to compute the bounds Ukip) and Lk(p) for the same k but 
many different values of p, the Hj^s need to be computed just once, and then 
they can be saved and used for all the different values of p. So, using (1.23) 
to compute Uk{p) for each of the many values of p, one will need to perform 
just 0{k) operations — raising p, . . . ,p + k to the powers -ffo,fc, • ■ ■ ,-fffc,fc or, 
what is essentially the same, compute 0{k) logarithms, namely, hip, . . . , ln(p + 
k) — followed by 0{k) multiplications and additions. In some computational 
environments (such as Mathematica) , the calculation of the logarithm of a given 
number may be about as fast as that the addition of two real numbers. Thus, one 
will need just 0(log j) (rather than 0(log^ i)) operations to make the relative 
errors (5io,fc(p) and (5up,fc(p) less than e £ (0, 1). 

There are a very large number of upper and lower bounds on the Wallis ratio 
and on the more general ratios of the form r(a; + s)/T{x), as well as related 
bounds for the digamma function tp, as in (1.15). The survey by Qi [13] lists 
204 sources (only a rather small minority of which concern just background 
information). See e.g. [2, 3, 5, 6, 7, 10, 16, 18]. However, there appear to be 
very few known convergent sequences of upper and lower bounds such as those 
presented in Theorem 1.1 and Proposition 1.3 here. 

A convergent sequence of lower bounds on the Wallis ratio (and, hence, by 
the identity (1.11), a convergent sequence of upper bounds) can be obtained 



immediately from the identity W{x) ~ y 2-F'i(— i, — a;; 1) a; due to Watson 
[18], which is in turn a corollary of a theorem by Gauss [1, Theorem 2.2.2] 
for the hypergeometric function. However, the hypergeometric series converges 
rather slowly, so that the number of operations it takes to make the relative 
errors less than e <£ (0, 1) using the Gauss-Watson identity is on the order of 
^2.^i/(2;+2)^ which, for any fixed x > and £ J, 0, is asymptotically much 

greater than the mentioned numbers log^ - and log -. Also, with our approach 
one has the extra bonus of the upper and lower bounds on the difference 
ip{x + 1) — ip{x + ^), as provided by Proposition 1.3. 

One should also mention the sequences of lower and upper bounds ak{x, s) := 
{x + s)'^-'{x + s + k- l)('=V(a; + fc)^''^ and /3fc(x,s) := {x + k + s^-^ix + s + 
k — V)^^'^ l(x + A:)('°) on the more general ratio p|||i[|]^j for any s G (0, 1) given 
in [15] (the same bounds, in different notation, were given also in [9]); however, 
_ 1 = (2i|±^) - 1 ~ as fc ^ cx), so that at least one of the two 

relative errors when using this result is on the order of \ and hence the required 
number of operations is very large, on the order of 7. 

In [16, 8], asymptotic expansions for W{x) are given, with the relative errors 
diverging in fc, but converging to as a; 00. Also, the sequences of bounds 
given in [14] are shown in [15] to diverge monotonically in k away from the 
target. 

Sequences of upper and lower bounds on t\}(x) were given in [9, (36)], where 
the difference between the /cth upper bound and the fcth lower bound is on the 
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order of 



, which is much greater tlian the geometricaUy decreasing absolute 



errors of the upper and lower bounds Uk{x) and £k{x) on il}{x + 1) — ip{x + i) 
in (1.18). 

Using methods similar to the ones presented here, one can obtain sequences 
of upper and lower bounds on the more general expressions of the forms r{x + 
s)/T{x) and tjj{x + s) — i^{x) for s € (0, 1), which converge faster to the target 
than the ones in the existing literature briefly reviewed here, but not as fast as 
geometrically. These results will be given elsewhere. 

Graphs of the (signed) relative errors '^''^^^7^^^^ ^'^'^ '^''''^(^^^'^ (*"^- (l-^^) 
and (1.14)) for k G {1, 2, 3} and x G (—5, 16] are shown as the matrix of plots 
in Figure 1; the columns of this matrix correspond to the values of k (shown in 
the plot labels), while the rows correspond to the different intervals of values of 
x: (-i,2], (2,6], and (6,16]. 




-0.0002 

-0.0004 










Fig 1. The relative errors of the hounds Ck{x) and Uk(x) for k G {1, 2, 3} and 2: S {— 4, 16] . 
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It is seen that the relative errors deerease fast in k and in x] moreover, they 
are rather decent even for fc = 1 and x £ (— i, 2]. Note also that, in view of the 
identity 



W{x) = W{x + m) 



jx + 



{x + 

for all natural m and all a; e (— ^, oo), without loss of generality one may assume 
X to be large. 

2. Proofs 

Here we shall prove Proposition 1.4, Proposition 1.3, and Theorem 1.1 - in this 
order. 

Proof of Proposition 1.4- Using the identity 

C^C;:.')^ ('";') 

(with the usual conventions that (J^Si^) = ("m^) ~ '-'^^ '^^^ ^^^^ Cp^m = 
Cp,m-i ~ fp+i_m_i for all natural to, where crp_,„ is the same as in (1.5), with 
Cp,o = hip. In turn, this leads to the identity 

Sp.k = I In ^ + 5 {Sp,k-i - Sp+i^k-i) (2.2) 

for all natural k, where Sp^k '■= lnJ7fc(p) = X]m=i '""'p,™' with Sp^ = 0. 
The identity (2.2) can be rewritten as 



""'^'^^i v^iuk-.ip+iy 

which is obviously equivalent to the first identity in (1.21). Also, in view of (1.6), 
(2.3) can be rewritten as the first identity in (1.20). Using the latter identity 
withp+1 instead ofp to replace Uk{p+1) in (1.6) by \/Uk-i{p + l)Lk-i{p + 1), 
one obtains the second identity in (1.20). In the latter identity, use again (1.6) 

to replace Uk-i{p + 1) by ^j^/Lk-i{p); this, together with the definition of 

Mk{p) in (1.21), leads to the second identity in (1.21). To check the "initial" 
conditions is straightforward. Thus, Proposition 1.4 is proved. □ 

Proof of Proposition 1.3. The four inequalities in (1.16) trivially follow from 
the definitions (1.17). Relations (1.18) follow by (1.16), (1.1), and (1.12), which 
imply il}{x + l) — ^{x + \) = i + 2-^^-1^^^. Thus, to complete the proof of Propo- 
sition 1.3, it suffices to verify the inequalities (1.19) (which will, in particular, 
imply the equalities in (1.16)). 
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Using the Dirichlet formula [1, Theorem 1.6.1] and then the substitution 
V := jiji 0116 has 

for any positive x and 5. Now the substitution t := v^^^ followed by an integra- 
tion by parts yields 



dlnr(p) _ r /p+l\ , fP\] _ ^ pf"^^ 



2p =p 

dp 



-1 = 2/ "^^^ - 1 = 2Jp,o, (2.5) 



where 



Here and subsequently, it is assumed that p G (0, oo) and k e {0, 1, . . . }, unless 
otherwise indicated. Further integrating by parts, one obtains the recurrence 
relation 

_ + (2 + fc) Jp.fc+i 

p + k + l 

whence, by induction, 

Jpfl =pik{p) +pRp,k, 
where £k{p) is as in (1.17) and 

(fc + 1)! 

^f^'^ p(fc+i) "'P'*^" ^^-"^ 

So, by (2.5), one has 

^^-4(p) = i?M. (2.8) 

For any t G (0, 1) and k ^ 0, the integrand i^^*^^yi-^k in (2.6) does not exceed 
(tti)'', which is less than 2^''. Hence, < i?„ i. < 2^^ — > and hence, 
by (2.8), 4(p) — ^ dXnrip) ^ ^j^g^^ j^g^g ^j^^ g^.g^ equality in (1.16). On the 
other hand, (1.10) implies 



1 dlnr(p+ 1) _ dlnr(p) 



2p{p +1) dp dp 

Therefore and because £k{p + 1) ^ d\nr{p+i) ^ definition of Uk{p) in (1.17) 
yields Uk{p) — > ^ 'dp''^'' ' as A; oo. In turn, this implies the second equality in 
(1.16). The first inequality in (1.19) follows immediately from (1.16). 

Next, observe that, by (2.6), Jp^k is decreasing in p. Since the denominator 
of the fraction in (2.7) is positive and increasing in p > 0, one sees that the 
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"remainder" Rp^k is decreasing in p > 0. So, in view of the definition of Uk{p) 
in (1.17) and the identities (2.9) and (2.8), 

, , dlnr(n) dhir(w+l) „ , , „ „ dhir(n) „ , , 

uk{p) ^ --lk{p+l) - Rp+i,k < Rp,k = — T^-h{p), 

dp dp dp 

whence the second inequaUty in (1.19) follows. 

The last inequality in (1.19) is obvious. Finally, introduce am ■= ^^(^+1™' 
and note that < i for all natural m, so that 

dlnr(p) _ ^^^^^ ^ ioo{p) - h{p) = ^ a„ < 2aA;+i, 

^ ni=k+l 

which proves the third, penultimate inequality in (1.19) and thus completes the 
proof of Proposition 1.3. □ 

Proof of Theorem 1.1. Take any p S (0, oo) and k G {0, 1, . . . }. Observe that 

= X^JLo ' — p-f j ' ^ foi' natural m; this can be checked by induction in 
m or as follows: 



noo ni 



J=0 

e-P"(l-e~")'"dM 

r(p)r(77i + 1) _ m! 



[\p-\i-zY 

Jo 



dz 



r{p + m + 1) ■ 



So, by (1.17), 



m=l j=0 ' ■' 



Recalling now (1.4) and (1.5), one has 

4(s)ds = -lnC/fc(p); (2.10) 



here, it was taken into account that X]JLo(~-'^)"' ('") ^ *-* ^'^"^ hence 
Er=o("lP(7) + j) = E;"o(-1P(7) Ms +3) - ln.s] 0, for each 

m = 1, 2, .... In view of (1.17) and (1.6), one similarly has 

Uk{s)ds^ TT^—TT- 4(S + I)ds 

Jp Ms + 1 ) Jp ^2 11) 

= - hi J^-r + hi Uk{p + 1) = - hiLfe(p). 
V 1 
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Also, it is easy to see that r{p) — )■ 1 &s p ^ oo. Therefore, 

ln.(p) = _ r^Mflds. (2.12) 
Jp ds 

So, the inequalities (1.7) follow by integrating those in (1.16) in view of (2.10), 
(2.11), and (2.12). Similarly integrating inequalities (1.19), one sees that 

0<ln4:(4<ln^<(^ + l)' r ds 



Lk{p) r{p) 2^+1 J J, s('^+2) 

2-'^-i(fc + l)! p ds 



<- 



(^ + 2)^=) Jp s(s + l) 

. 2-"-i(fc + l)! , P+1 , 
~ (p + 2)W ^^^<^f^^' 

where p* is as in (1.9). This yields (1.8) and also implies that Lk{p) and Uk{p) 
both converge to r{p) as fc — > cx), so that the identities (1.3) follow as well. Thus, 
Theorem 1.1 is proved. □ 
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